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Abstract 

Potential algebras can be used effectively in the analysis of quantum systems described by (2 + 
1) dimensional Dirac equation. We find integrals of motion of Dirac operator which close centrally 
extended so(3), so(2, 1) or oscillator algebra. The algebraic framework is used in construction of 
physically interesting solvable models. It is applied in description of open-cage fullerenes where the 
energies and wave functions of low-energy charge carriers are computed. Besides, the potential algebras 
are employed in the context of one-dimensional, both relativistic and nonrelativistic, shape-invariant, 
quantum systems. We show that shape- invariance of the considered Dirac operators is associated with 
the nonlinear N = 2 supersymmetry. The potential algebra is used in construction of nonrelativistic 
shape-invariant systems with position dependent mass. We demonstrate that solvability of these models 



O j is associated with their N = 8 supersymmetry. 

1 Introduction 

Exactly solvable models play exceptional role in physics. They are simple enough to be solvable, yet 
sufficiently complicated to grab the essence of physical reality. In quantum mechanics, exact solvability 
usually means that the eigenfunctions and eigenvalues of the Hamiltonian can be found explicitly. They 
can be computed analytically by solving directly corresponding differential equations. 

Alternatively, they can be obtained in algebraic manner with the use of the integrals of motion; 
when the symmetries of the quantum system form Lie algebra, the spectrum and eigenstates of the 
energy operator can be deduced from the analysis of the admissible unitary representations [T], [2]. The 
symmetry operators preserve domain of the Hamiltonian by transforming one physical state into another. 
When they are time-dependent, they do not commute with the energy operator and, hence, they do 
not preserve energy. The associated algebra of the integrals of motion is called dynamical or spectrum 
generating, see [3], [3]. When the integrals of motion commute with the Hamiltonian, they rather reflect 
spectral degeneracy of the system. The algebraic structure, where the Hamiltonian plays the role of 
central element, is denoted as potential algebra. In the current article, we will focus on the analysis and 
applications of the potential algebras of the relativistic quantum systems described by (2 + 1) dimensional 
Dirac equation. 

The low-dimensional Dirac equation appears in surprising variety of physical settings that are studied 
both in quantum field theory and in the condensed matter physics. Let us mention the (1 + 1) dimensional 
variant of the famous Nambu-Jona-Lasinio (chiral Gross-Neveu) model [5], [5], [7J, analysis of fractionally 
charged solitons [8], [9] as well as description of inhomogeneous superconductors [10], carbon nanotubes 
|11| . |12j . or linear molecules like polyacetylene [13] . The (2 + 1) dimensional massless Dirac Hamiltonian 
appears in the low-energy approximation of dynamics of charge carriers in graphene and in related carbon 
nanostructures [2], [15], [16], [T7J. 

The work is organized as follows: in the next section, the potential algebras are constructed for a 
wide class of Dirac operators with the use of a specific ansatz for both the ladder operators and for 
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the structure of the algebra. Relevant aspects of the representations based on the lowest/highest weight 
vectors are discussed in more detail, as well as the application of the algebraic structure in description 
of Dirac fermions in the curved space. The results are used directly in the third section, where solvable 
models describing low-energy charge carriers on the surface of open-cage fullerenes are analyzed. Two 
configurations are considered, with one and two holes in the surface of the crystal. In the fourth section, 
we apply the potential algebras in construction of one-dimensional shape-invariant systems. First, we 
focus on the shape-invariance of one-dimensional Dirac Hamiltonians and show that it is associated 
with N = 2 nonlinear supersymmetry. Next, nonrelativistic systems with position dependent mass are 
studied. The algebraic background of the shape-invariant superextended Hamiltonian is identified as 
N = 8 superalgebra. Illustrative example of solvable system is presented. In the last section, we discuss 
briefly two specific two-dimensional systems; we make few comments on Dirac oscillator and its dynamical 
symmetries. Additionally, we show that the so(2, 1) potential algebra can be used for analysis of (1 + 1)- 
dimensional quantum system whose metric can be identified with the (restricted) BTZ black hole metric. 
We conclude with short discussion of the results and outlook to possible future research. 



2 Potential algebras of two-dimensional Dirac operator 

Let us consider the operator ho which is given in terms of the Pauli matrices a a and real functions g\, 
g 2 and g 3 , 

hD = io\ (gid Xl + + (T 2 {-ig2d X2 + C353) , g a = g a (xi), a = 1,2, 3. (1) 

We leave the domain oihn as well as the range of the real variables x\ and X2 unspecified at the moment. If 
not stated otherwise, the coordinates x\ and x 2 are considered to be space-like. Then ([T]) can be regarded 
as Dirac Hamiltonian of (2 + l)-dimensional Dirac equation. The function g 3 = g%{x\) represents external 
magnetic field with real coupling constant C3. The operator is manifestly hermitian with respect to the 
standard scalar product. It anticommutes with 03, which implies that its spectrum will be symmetric 
with respect to zero. 

The system represented by ([T]) has integral of motion J3, 

J 3 = -id X2 , [h D ,J 3 }=0. (2) 

We are interested in the settings where two additional integrals of motion J± exist, such that 

[J 3 ,J±]~±J±, [J_,J + ]~J 3 , [J±,h D }=0. (3) 

Let us make the ansatz in the form J + = e tx ' 2 fX)fc=i(Afe,l+ a 3 J ^k,2)d Xk +C\ + 03C2) where the coeffi- 
cients ^4fc,a an d C a (a, k = 1, 2) depend on x\ only. The operator J_ is defined as J_ = j\. Inserting the 
ansatz into ([3|) and comparing the coefficients of corresponding derivatives, we find that the coefficients 
of J± can be expressed in terms of the functions g±, g 2 and g 3 as 

t ■ ±1x2 (a 1 9x 1 ffiff 2 (± J 3 + |) gig' 3 a 2 \ ,,x 

J± = ie^ gi d Xl + — H — ± =F 170-3 • (4) 

\ 2 92 92 2 J 

It suggests that the function g 2 should be node-less to avoid singularities in definition of these operators. 
Additionally, the functions g%, g 2 and g 3 have to solve the following set of differential equations, 

9293 + 9\929z ~ 9i92(g'i93 + 5103 ) = °> 

4,2/2 /II, " \ n 

92 +9192 - 9192(9192 +9192) = °> 
(91 {92 (9192)' ~ 2ffi9 2 2 )\ 



9l 



0. 



' 91 (92(9193)' - ZgigWz) V 
g\ J 
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Quick inspection suggests that the system can be reduced to two equations as long as 53 = 92- However, 
we will omit this solution; looking at jl]), it would just shift the operator J3 by a constant C3. The solution 
of the unreduced system can be found such that g\ as well as 53 are given in terms of 52 in the following 
manner, 

n A , r n 2 [ 4(92 +ci), for ci ^0, 

9l = 9l = (5) 

072 } ( c\g- 2 \ for Cl = 0, 

where the c\ and C3 have to be real valued to maintain hermicity of the Hamiltonian. 

The operators J3 and J± form Lie algebra. Its explicit form depends on the value of c±. For c\ 7^ 0, 
it can be identified either as so(3) (ci < 0) or so(2, 1) (c\ > 0), while the oscillator algebra is restored for 
ci = 0. ho commutes with both J3 and J±, and, hence, it is the central element of the following algebra 

[^,J ± ] = [^,J 3 ] = 0. (7) 

Let us notice that nonrelativistic systems with similar algebraic background were classified as superin- 
tegrable in [18J. The representation space of ©-([71) is spanned by the eigenvectors of ho- It can be 
decomposed into subspaces where ho acquires constant value and where the algebra ©-(0) is irreducible. 
These subspaces can have either finite or infinite dimension, dependently on the type of the algebra. 

The square of ([T]) can be written as a second-order polynomial in the generators of ©. When c% 7^ 0, 
it coincides up to an additive constant with the standard Casimir operator C = J+ J_ — c\J^{J^ — 1) = 
J_ J + — c\J^{J-i + 1) of the algebra so(3) or so(2, 1). When c\ = 0, the algebra is not semi-simple (Killing 
form is degenerate) and the standard (quadratic) Casimir operator does not exist. There holds 

^ = f J + J_-c 1 ((j 3 -i) 2 -c|)=J_J + -c 1 ((J3 + i) 2 -c|) for Cl /0, 
D \ J+ J_ + 4c 3 (J 3 - \) = J-J++ 4c 3 (J 3 + |) for d = 0. 

The operator commutes withQ 1Z = R X2 R cs a\. The reflection operators R X2 and R cs are defined as 
Rx 2 f(xi, X2, C3) = f(xi,— X2,cs) and R C3 f(xi, X2, c 3 ) = /(xi, X2, — c 3 ). For generators of the potential 
algebra, there holds 

izj-iz = j+, tzj 3 tz = -j 3 . 

In our analysis, we will focus on the representations that can be constructed from the lowest (highest) 
weight vectors, the zero-modes of ladder operator J_ (or J + ). These vectors, denoted as (or V'm); can 
be fixed as the mutual eigenstates of hr> and J3, 

J±^m = 0, J^t = m ^m, h D 1p± = E±lp±. (9) 

The energies can be found with the use of the relations (JSJ , 



{ Eif = { -ci((™±h) 2 -4) for Cl /0, (J()) 
4c3 (m ± 2) for ci = 0. 



To keep the energies real, there must hold \m ± | > C3 for ci < 0. For ci > 0, we get mrb M < C3. 

The other eigenstates of hp corresponding to the same energy can be found by repeated action 
of the ladder operator J + (J_) on the vectors Vm ( or ^m)- These eigenstates then establish irreducible 
representations of ©-([7J that is specified by the energies of the corresponding lowest (highest) weight 



1 Similar operator that switches the sing of coupling constant was employed in [20] in the analysis of quantum systems 
with nonlinear supersymmetry. 
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vector. These representations are usually associated with the bound states of the quantum system, see 
e.g. P2|. 

The unitary representations of so(3) are finite dimensional. Consequently, the lowest weight vector, 
let us fix it as ipZ m , has be to annihilated by specific power of the ladder operator J + . There holds 
(J+) d ipZ m = where d denotes dimension of the corresponding representation. The quantum number m 
has to be either integer or semi- integer; we have ip^ ~ ( J + ) d ~ lr il>Z m so that J^JjI~ lr 4>Z m = mJ^V-m' 
Additionally, we can use © and write J3 J d ~ 1 i^Z m = (— m + d — 1) J+ _1 VC m - Combining these relations 
together, there must hold —2m + d — 1 = 0. Since d is a positive integer, we get that m has to be 
non-negative integer or semi-integer. The admissible values of the coupling constant C3 depend on m; 
when m acquires integer (half-integer) values, C3 has to be half- integer (integer). We refer to Appendix 
A for details of the proof. Let us notice in this context that the representations of odd dimension for two 
different half-integer values of C3 were considered in j2TJ. Representations of so(3) and so(2, 1) were also 
discussed in [22] where they were obtained analytically with the use of the master equation. 

In the end of the section, we shall discuss shortly the situation when the spin— 1/2 particle is living in 
the two-dimensional curved space. Let us consider the space where the metric tensor acquires diagonal 
form and the non-vanishing components depend on one coordinate only, 

„ _ ( #uOi) 

9vv ~ \ g 22 ( Xl ) 

In the most of the text, we will deal with the systems where both x% and x 2 are considered to be 
space-like (gn > 0, 522 > 0), i.e. the metric tensor corresponds to a curved surface embedded into 
three-dimensional euclidean space. Quantum settings of this kind appear frequently in the analysis of 
carbon nanostructures where deformations of graphene crystal are considered |24| . The massless Dirac 
Hamiltonian can be written in the following form (see Appendix B for details), 

hn = ^l-TTZ^-TTZ + a i—T= + C39302- (11) 

It is manifestly hermitian with respect to the standard scalar product. The operator (|lip is a special case 
of (H|); when g\ and g 2 in ([!]) are both positive, we can fix g\ = -7=, 52 = T/j^y- The operators ([TJ and 
(|lip then coincide. Employing the same substitution in and (0), we find that the Hamiltonian (jlip 
has the integrals of motion 

provided that 511 and the external potential 53 are fixed in the following way 



a ~ 1 (g ^ 2 „ 

5n — 7 77— r, 9s 

4^22(1 + cig 22 ) 



C3 522 2 ' for ci = 0. 



(13) 



Concluding the section, let us notice that the generally non-constant coefficients of d xi and d X2 in (jlip can 
be interpreted as the space-dependent modulation of Fermi velocity of massless Dirac particle in graphene 
induced by the strains and ripples of the crystal [25J. 



3 Solvable models of open-cage fullerenes 

Fullerenes are spherical molecules made of carbon atoms. The spherical shape is due to twelve pentagons 
that are inserted into hexagonal carbon lattice. The most famous example is the Cqq, "Buckminster- 
fullerene" [26J, which can be classified as a truncated icosahedron with 60 vortices and 32 faces. Fixing 
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the number of pentagons in the lattice while increasing the number of hexagons, we can get not only 
bigger fullerenes like C240 or C540, but also objects of more complicated shapes, e.g. capped nanotubes, 
elliptic fullerenes etc.. However, we will be interested just in the systems where the surface of the crystal 
can be approximated by the spherical geometry. 

In quantum chemistry, the situation is considered where small enough atoms or molecules (e.g. of 
hydrogen) are inserted inside the fullerenes. During the process called "chemical surgery", some bonds 
between the carbon atoms of the fullerene are broken, making the hole into the surface. The alien atoms 
are then inserted into the opened cage and the hole is closed again [27]. In this section, we will focus 
on the electronic properties of the open-cage fullerenes with one or two holes in the surface. We will 
construct exactly solvable models that can serve as the low-energy approximation of these systems. 



3.1 Spectrum of fullerenes 

First, let us go briefly through the analysis of the electronic properties of the (unopened) fullerenes. The 
low-energy excitations of charge carriers in the crystal can be investigated efficiently within the framework 
of the Dirac equation on the spherical surface [21]. The surface embedded into three-dimensional space 
can be parametrized in spherical coordinates as x = sinxi COSX2, V = sinxisinx2, z = cosrci where 
x\ E (0, 7r) and X2 = (0,2ir). The associated metric tensor is explicitly 

( 1 

9txv ~ \ sin 2 xi 

The entries of the metric tensor satisfy (|13p for c\ = — 1. Employing the formula (jlip . we can write down 
immediately the Dirac Hamiltonian corresponding to the spin— 1/2 particle on the sphere, 

hD = i&id Xl + fJ 2 - 1- C3<7 2 cotxi, xi G (0, ir). (14) 

smxi 

The potential term corresponds to Dirac monopole that is situated in the center of the sphere. It occurs 
due to the topological defects, twelve pentagons, that have to be inserted into the hexagonal lattice of 
graphene to close the spherical surface [16] , [T7J , [21] . The actual value of the coupling constant is C3 = | 
where n corresponds to the number of conical defects in the crystal. As n = 12 for fullerenes, we get 
c 3 = §, see US], |2U. 

We will not write down the explicit form of J±, it can be extracted easily from (|12p . The lowest weight 
vectors can be found explicitly by solving Q. They are 
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, - e*- (sin*!)-** ( ) , WF = 1^1, , ,.5) 

h D^ = \l [ m - 7; ) - 4 (16) 

We require the eigenfunctions of (|14p to be square integrable and vanishing at x\ = 0, tt. The vectors (|15p 
are physically acceptable as long as m < — \c% \ + \ . In that case, they are regular and correspond to real 
energies, see the restriction below (|10p . Degeneracy of the zero-energy level is equal to the dimension of 
the representation of so(3) spanned by zero modes of the Hamiltonian. For C3 = |, the zero-energy has 
triple degeneracy. 

For purposes of the forthcoming analysis, it is convenient to define transformation G F 

' 1/2 x 



G F = tup 1 ^sinxi [ f _ 1/2 ] 01 
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and rewrite the lowest weight vector in the factorized form 



^-=e^y- m G F \ V ? y = mn Xl , w F (y) = { (17) 

where 

k u = 1 - 2c 3 , fc d = -2c 3 - 1. (18) 
With the use of this notation, the ladder operator can be rewritten as 

G- l J+G = ie iX2 (d Xl -cot Xl J 3 )= ie iX2 (l-yw F (y))(d y -^-Y (19) 

It would be necessary to check that a multiple action of J+ on ^~ does not violate regularity of the wave 
functions. However, we can skip this task as the considered wave functions are known to be regular and 
based on the Jacobi polynomials, see [22], [25]. It means that the action of (fTUj) on (fTTj) preserves the 
required boundary conditions. We will use this fact in the analysis of the open-cage fullerenes. 

3.2 Open-cage fullerene with one hole 

Let us consider the fullerene molecule with one hole carved into its surface (see Fig.l). We suppose that 
the twelve pentagons allowing the crystal to close the spherical shell are still present in the lattice. To 
meet this requirement, we can suppose that the model describes larger molecules like C240 where one or 
few hexagons can be extracted from the crystal without affecting the pentagons. The considered spherical 
surface can be parametrized as 

x = sm(irdn(xi,k)) COSX2, y = sin (irdn(xi, k)) sinx2 z = cos (irdn(xi, k)) , 

where x\ £ (0,K), x<i = (— 7r, it). The components of the associated metric tensor 



(nk 2 cn(xi, k)sn(xi, A;)) 2 

sin 2 (ftdn^xi, k)) 



satisfy (|T3|) for c% = — 1. 

The functions dn(xi,k), cn{x\,k) and sn(xi,k) are Jacobi elliptic functions and K = f Q 2 (1 — 
fe 2 sin 2 6) _1 dO is the complete elliptic integral of the first kind with < k < 1. Let us briefly refresh rele- 
vant properties of these functions. dn(x\, k) is periodic with periodicity 2K while cn(x\, k) and sn{x\, k) 
have periodicity 4K. There holds < dn(xi, k) < 1 for all real x\ and all values of the modular parameter 
k £ (0,1). The functions cn(xi,k) and sn(xi,k) can be understood as a generalization of trigonomet- 
ric functions; they satisfy the relation cn(xi,k) 2 + sn(xi,k) 2 = 1. There holds cn(xi,0) = cos 21 and 
sn(xi, 0) = sinxi, dn(xi, 0) = 1. When k — > 1, the period of the Jacobi elliptic functions tends to infinity. 
There holds sn{x\, 1) = tanhrci, whereas sn{x\, 1) = dn{x\, 1) = sechxi. In the cases where the actual 
value of the modular parameter is not informative or it is clear from the context, we will use the shortened 
notation sn(xi) = sn(xi,k), cn(xi) = cn(xi,k) and dn{x\) = dn(xi,k). 

Dirac Hamiltonian on the considered surface acquires the following explicit form 

1 / dnxi{sn 2 x\ — cn 2 xi)\ J3 

h D = io\—, d Xl H +cr 2 . r + c 3 (T2Cot(7rdnxi), (20) 

k z TTcnx\snxi \ lcnx\snx\ J sm(7ran:ci) 

where x\ G (0,K). Likewise in (|14p . the potential term corresponds to the field generated by Dirac 
monopole. We will require the eigenstates of ho to be regular and vanishing at the edge of the hole (i.e. 
at x\ = K). The lowest weight vectors f/'m annihilated by J_ are 



1/2 



w l \ A/ c 3 - ±+mWj 
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where WI = and h D ip m = ^{m - 1/2) 2 - c^ m . 

The wave functions (|2ip comply with the required boundary condition (cn{K) = 0). However, it is 
not clear whether the multiple action of J+ on (|21|) does not violate the boundary conditions. We can 
analyze this point using the close analogy with the model of fullerenes. Let us define the transformation 
Gj as 

' -1/2 ^ 



Gi = \/ cnx\snx\ sin(7rcZn(xi)) Wj 1 1 1 / 2 

y o Wj 

The lowest weight vector can be then written as 



« = <~V~G, I V 5t±i^"f , „ _ sin(mlnXl) , „,(„) _ J ^" 2 > (22) 

'ca-^ + mwJ* J { + y > dnxi>l/2, 

where the constants k u and kj are those in (|18p . Performing the "gauge" transformation of J+, we get 

Gfj+Gj = ie iX2 (- 2 8 X1 + cot Trdn Xl J 3 ) = -ie lX2 (l -y Wl {y)) ( d y - —) . (23) 

yk^ircnxisnxi J \ y J 

The latter operator is identical with (fP9"j) . The vector (f2"2"j) coincides with (fT7|) up to gauge factors Gf and 
Gi and the up to the sign of the upper component. Hence, likewise in the case of fullerenes, the multiple 
action of J + on ^~ will keep regularity of the wave functions, i.e. only positive powers of y will emerge. 

3.3 Open-cage fullerene with two holes 

Similar results can be obtained for spin-1/2 particle on the spherical surface with two holes centered in 
the poles, see Fig.l. The surface can be parametrized as x = dnx\ cosx2, y = dnx\ sinx2, z = ksnx\, 
where x\ G {—K,K) and X2 £ (—tt,tt). The metric tensor reads this time 

/ (kcnxi) 2 
9 ^-\ dn 2 Xl 

The explicit form of the Hamiltonian can be written as follows 

iai ( dnx x snx x \ J 3 ksnxi 

h D = -, I d xi H +^23 \-c3a2-; , x 1 £(-K,K). (24) 

kcnxi \ lcnx\ J dnx\ dnx\ 



Figure 1: The models of open cage fullerenes. We fixed k = .99 in the left figure while k = .9 in the right 
figure. 
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Its eigenvectors should be regular and vanishing at the edges of the holes (i.e. for x\ —> ±K). The lowest 
weight vector which satisfies ho^m = y (m — \) 2 — nas the explicit form 



^ = +2 l i — r— -.,,-1/2 > w " = H^— )• ( 25 ) 



W 



ii ^ Wc 3 -2+m%' 

It complies with the required boundary conditions since cn(K) = cn(—K) = 0. In order to show that 
application of J + on these vectors does not violate the required properties, let us define the matrix 



Gjj = yj cn{xi)dn{x\) vffj ' " " 



W TT 



II n 1/2 

1 Wjj 



The lowest weight vector can be written in the following form 



= e^y- m G n \ \ ), y = dn Xl , w n {y) = <j i+ J^' \ ' (26) 



c 3 - ± + m Wj| 
while the ladder operator is transformed into 

GjtJ+Gu = uf°* ( d Xl + = _ ie - 2 (i _ ywil{y) ) U _ J A . (27) 

\kcnx\ dnx\ J \ y J 

Comparing ([26]) with ([T7]) and ([27]) with (fT9p . we can conclude that the repeated action of J+ on ip^ will 
preserve regularity of the wave functions. 

The formulas ()22|) . (|23p . (|26p and ()27p allow to find the required eigenstates of the corresponding 
Hamiltonian as well as energies which are given by (jlOp for c\ = — 1. Each of the systems has distinct 
topological genus as the number of the holes in the spherical surface varies. Despite this distinct topological 
nature, the models are spectrally identical due to the same algebraic background. 

In the next section, we will discuss a different application of the potential algebras. We will con- 
sider one-dimensional, both relativistic and non- relativistic supersymmetric systems that possess shape 
invar iance. 



4 Shape invariant Dirac Hamiltonians and PDM Schrodinger opera- 
tors 

Shape invariance underlies exact solvability of the most of the one- dimensional nonrelativistic quantum 
mechanical systems. It originates in supersymmetric quantum mechanics based on the Witten model [29], 
[30] . [31] , Within this framework, two hermitian Hamiltonians, 

Ht = "C + Vm = -C + W(xx, mf ± W'( X1 , m), (28) 

are intertwined by the operator q m = —id xi — iW(xi,m) and its conjugate q}n, i.e. H^qm = q m Hm and 
H~ qti = qla H+. The function W(x\,m) is called superpotential. In the context of the analysis of differ- 
ential equations, the intertwining operators are known as Darboux transformations [32]. Supersymmetric 
system is defined in terms of the extended matrix Hamiltonian H m and integrals of motion Qm and Qm , 

Hm = -d 2 Xl +W(x 1 ,m) 2 + a 3 W'(x 1 ,m), Q$ = -iaid xi +a 2 W(x 1 ,m), Q$=mx 3 Q«. (29) 

In close analogy to the supersymmetry in quantum field theory, these operators compose a superalgebra. 
The generators (|29p are classified as bosonic or fermionic provided that they commute or anticommute 
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with a grading operator. In our case, the grading operator can be identified as <j 3 . Then H m is bosonic 
([H m ,as] = 0) and Qm are fermionic operators ({Qm ^3} = 0, a = 1,2). The Hamiltonian H m then 
possesses N = 2 supersymmetry associated with the superalgebra 

= 0, {Q%\Q<$} = 28 ab H m . (30) 

Let us suppose that the potentials and V~ differ just in the value of the coupling constant and in 
additional constant term, i.e. V£ = V m +$ + r(m) where 5 and r{m) are real. Then the Hamiltonian 
is called shape-invariant as V,^ and are qualitatively of the same shape. This property can be used 
to compute spectrum and wave functions of very effectively just from the knowledge of the ground 
state, see e.g. [30]. We will discuss this point in more detailed in the section B~2l 

The condition of shape-invariance (|3ip can be written in terms of the function W(xi,m) in the 
following manner 

W{x 1} mf + W(x 1: m)' = W(xi,m + 5) 2 - W(xi,m + 5)' + r(m). (31) 

This formula proved to be very useful for the systematic search of the shape-invariant systems which has 
been done in the literature extensively, see e.g. [33], [33] . 

Analysis of shape-invariant Dirac operators in the context of supersymmetry is more complicated. 
The framework based on the superpotential similar to (|28p - (|3ip is absent. Usually, shape-invariance of 
Dirac Hamiltonian is treated such that the square of the operator is identified with a known, nonrela- 
tivistic, shape-invariant operator after a series of transformations [35]. In this way, exact eigenstates and 
eigenvalues of the original relativistic system can be found. 

We will present an alternative approach inspired by Balantekin and Gangopadhyaya [36], [37]. They 
found that the shape-invariance of the nonrelativistic systems can be understood in terms of a higher 
dimensional Hamiltonian that possesses Lie algebra of integrals of motion. In this framework, the super- 
partner Hamiltonians as well as the associated intertwining operators q m and q m correspond to an 
appropriate restriction of the higher-dimensional Hamiltonian and of the ladder operators to the subspaces 
with fixed value of angular momentum. 



■"raj Vm 



4.1 Shape-invariance of Dirac operators via N = 2 nonlinear supersymmetry 

Let us consider restriction of hn in (pQ) and of the ladder operators (jH) to the subspaces where J3 acquires 
fixed value. The operator ho reduces to the one-dimensional Dirac Hamiltonian 

h m = h D \ J:i=m = e- imX2 h D e imX2 = ia x ( 9l d Xl + y) + + g 3 a 2 . (32) 

The ladder operators J± exchange the subspaces where J3 = m and J3 = m + 1. Corresponding restricted 
operators j± >m read explicitly 

j m = e -^)-J + e— = iLd Xl + f + ^±il + M_0, 

\ 2 92 92 ^ J 

Jm = e—U.e^^ = i Ld Xl + ^ - 9 ^ {m + - M + 92 \ ( 

y 2 92 92 ^ J 

Keeping in mind that [ho,J±] = 0, we can write down the following intertwining relations between 
the operators h m and h m+ i, 

hmj m = J m ^m+1) ^m+lJm = Jm^m- (^4) 

These relations establish shape-invariance of the one- dimensional Dirac operator h m mediated by the 
operators j^. They also make the link with the supersymmetric quantum mechanics. The operator 



9 



(j m ) can be identified as the generalized matrix Darboux transformation. General properties of these 
operators were discussed in |46j in detail. In particular, it was proved that the product j m j^i or 3^n3m ^ 
proportional to the second-order polynomials in h m or h m+ i respectively, 



3m3m ~ h"m ^"mi 3m jm ~ ^mi (35) 

where £ m = (E+) 2 . 

It is convenient to treat the system in the framework of the extended matrix operators, 



h m+1 J ' " m \ m+1 J \ 2 J 2 ' 

J ™ = {% ) ' J ™ = i7 *Jm> ^ = a z ®l = diag(l, -1), (36) 

\ Jm U J 

where 1 denotes the two-by-two identity matrix. The operators act on the bispinors with the upper spinor 
from the subspace where J3 = m and with the lower one from the subspace where J3 = m + 1. 

(3) 

The matrix T3 can be identified as the grading operator of the superalgebra, classifying h m and J m 
as bosonic while J7m and J7m as fermionic operators. The N = 2 superalgebra then reads 

{ Ji a) , J®} = 25 ab (h m - £ m ) , a, 6 = 1, 2, (38) 

where e a b is completely antisymmetric in indices and £12 = 1. The shape-invariance of the relativistic 
Hamiltonian (|32j) is encoded in the first commutator of (J37J). The superalgebra is nonlinear, see [38], as 
the anticommutator of the supercharges is a quadratic polynomial in the (extended) Hamiltonian h m . 
Hence, the shape invariance of h m is associated with nonlinear N = 2 super symmetry. 

4.2 Shape-invariance of nonrelativistic systems with position-dependent mass and 
N = 8 supersymmetry 

Another interesting application of the potential algebra (J6j)-(ITJ) can be done in the realm of the nonrela- 
tivistic quantum mechanics. The Schrddinger operators with the position-dependent mass M(x\), 

' -d Xl -L= d Xl — -±-j- A +V(x 1 ), (39) 



M(x x )i/4 X1 ^M{^) Xl Af(xi)V4 

emerge in condensed matter physics, e.g. in description of the semiconductors or heterostructures 
|40j . These systems have attracted lots of attention in the literature, see |41j-|42] and references therein. 

Let us pick up Dirac Hamiltonian (jlip . The square of the operator in the subspace of fixed angular 
momentum J3 = m can be written as 

1 „ 1 „ 1 f m \ 2 1/m x ' 



H m = m d *i ~l= d xi "T7Z + ~F= + c 3^3 - °3^= —== + C353 • (40) 



We can see that the diagonal entries of the matrix operator are of the type (|39|) provided that we fix 
M(x\) = g\i{x\). By construction, the Hamiltonian H m together with the supercharges Qm\ 

Qm = h m = icri^j7d Xl -^j- + a 2 -^ + c 3 g 3 a 2 , Q$ = ia 3 Q$ , (41) 
close the superalgebra, 

H m ,Q$] = 0, {q#,Q#} = 2S ab H m , (42) 
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which is identical with ([HP]) . Hence, (gDD-dMD extend the Witten's N = 2 supersymmetry dSSD"® to the 
systems with position dependent mass. 

The spectral properties of H m can be easily derived from those of Q$ • The spectrum of the latter 
operator is symmetric with respect to zero. It implies that the eigenvalues of H m have double degeneracy 
except the zero energy which is non-degenerated. Degeneracy of the levels is reflected by the supercharges 
Q$ ; they interchange the states of the same non-zero energy while they annihilate the zero mode of the 
Hamiltonian. 

The Hamiltonian H m is solvable due to the shape-invariance mediated by the operators j± m , 

Hm+ljm = im-^mi HmJm = Jm^rn+l- (43) 

These relations can be used to find the eigenstates of H m . Let us suppose that the ground state(s) of H m 
can be identified with ip + = -i/j+j j 3 = m (and gt3^> , in case that the ground state energy is degenerate). The 
ground state energy is £ m = (E m ) 2 . Then the first excited state of H m can be computed as jm^m+i where 
V'm+i is the ground state of the Hamiltonian H m+ \. The second excited can be computed 
while the n-th excited state can be computed as j m j m+ i ■ ■ ■ 3rn+n-i' l J J m+n w ^ ere ^™+ n * s ^ g roun d state 
of H m+n . The intertwining relations (|43p ensure that these functions are indeed the eigenstates of H m 
corresponding to the energy £ m +i, £ m +2 and £ m +n, respectively. 

The shape-invariance of H m expressed by (f4~3j) and the supersymmetry (f4~2j) can be unified in the 
extended supersymmetric framework where both Q$ and J$ play role of the fermionic supercharges. Let 
us mention that supersymmetric structure of a nonrelativistic system with two distinct sets of supercharges 
was discussed e.g. in [33] and [44] . It is convenient to define the following extended operators 

«-=(*" *!.)• c -(f qL)' <44) 

We fix CJ3 <%> (J3 as the grading operator which anticommutes with both Q$ and J$ . We define four 
additional fermionic operators T^Qm^ and T\Jm^ where T3 was fixed in (|3"6]) and T\ = 1<8>0"3 = diag{pz, 03). 
Additionally, we introduce four bosonic operators Xm' b ^ = Jm ] Qm ', a,b = 1,2. We have eight fermionic 
operators which we denote in the following manner, 

n(a,i) = q(o) n( a > 2 ) = T „n( a ) ji^ 1 ) = 7O) j( a < 2 ) = T , j( a ) 

•^"in -<"m ' -"^m — o-^-m u m ' J m ' u m ' L u m ' 

The set of eight bosonic operators, 

7~Lm, <Jrn\ ^m' i T l ^3) 

can be reduced to seven since the operator Jj, is not independent. It is a linear combination of identity 
matrix and of the operator T3, see (|36p . Hence, the commutators involving this operator are omitted. We 
also omit most of the relations which are identically zero. First, we can write 

[H m , Qti b) ] = 0, {Qt c) , Qt d) } = ^ ab (Scd + (1 - 5cd)r 3 )n m , (45) 

[Hm, Jt b) ] = 0, {jfr C \ J£ 4) ] = habited + (1 - 8 cd )n)(H m ~ £m), (46) 

[n, q£ c >] = -2ie ab Q£> c \ [r 3 , jM] = -Ke ah J^ . (47) 

2 The zero mode tpo of Q m ', Qm^ipo = 0, can be written as ipo = 9 3 ' 1 



a 2 e 



/9TT 



constants. If the upper component of the spinor is square integrable, the lower cannot have this property and, hence, we 
have to put Q2 = 0. 
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The sub-superalgebra (|45|) with c = d = 1 reflects the Witten's-type supersymmetry (|42p whereas the 
shape-invariance ()43|) of iif m is encoded in the commutator of (|46p . The superalgebra is completed by the 
relations which mix the supercharges Qm and Jm i giving rise to the bosonic operators Xm' b \ 

{Ji a ' c \ Q { t d) } = 25 cl 5 dl xt b) + 25 c2 5 d2 e ar e bs X^ s \ (48) 
[X<* b \ Xj***] = -2i (e bd (8 ac n - ie ac ) + e ac (S bd T 3 - ie db )) H m (H m - £ m ), (49) 

[Xt b \ Q^ d) ] = -2i (e bc 5 dl jt 2) ~ 5 bc 5 d2 e ar J^) H m , (50) 

[X$n\j£' d) ] = -2t (e ac 5 dl Q^ - 5^d d2£br Q^) (H m - £ m ) , (51) 

[T U Xfr h) \ = -2ie bc xt c \ [73, X<£' b) ] = -2ie bc x£ b \ a, b, c, d,r,s = 1, 2. (52) 

In the computation of the commutators and anticommutators, we used [Jm\Qm] = 0, [Jm\jm] = 
-2ie ab T ?J (U m - £ m ), [Qm,Qm] = -2ie ab T{K m . The superalgebra (@H])-([H2]) is nonlinear again. Its 
linearity is restored in subspaces of fixed energy where Ti m acquires constant value. The relations ([18])- 
([52]) establish N = 8 supersymmetry of the extended Hamiltonian % m with position-dependent mass. 



4.3 Illustrative example 

Let us present an example of a solvable shape-invariant system with position dependent mass. Since the 
effective mass in ([39]) is required to be bounded and positive for all x±, we fix the coefficients in (fll~]) as 

/ \ 2 / \ 2 CTh ( X 1 ) , , 

522 = sn[xi) , g n = dn(x 1 ), g 3 = c 3 — ? — r-, x\ G (0,2it). 

sn(xi) 

It complies with (|13j) for ci = —1. Hence, we deal with quantum system which possesses so(3) potential 
algebra. The explicit form of the supercharges (|4ip reads 



1 n 1 w cn(xi) ,_„. 

V^(xi) \/dn(xi) sn(x{) sn(xi) 

The nonrelativistic Hamiltonian acquires the following form 

H = 1 d , 1 d 1 | / m + c 3 cn(xi) \ 2 | ^ c 3 + mcn(xi) 
dn(xi)5 Xl dn(xi) 8,1 ^(^Jg \ snfci) / 3 s^(^i) 2 

The effective mass of the nonrelativistic particle coincides with dn 2 (xi) and is strictly positiv^. The 
particle lives in the finite interval in presence of the external potential. We require the wave functions of 
H m to be regular and vanishing at the borders of the interval. 
The explicit form of the operators j m ± reads 

it = i { -rJ=d Xl ^J= T (m+ h ] m } Xl) t4,T ^h~\ ■ (55) 
\^M^) 1 ^/M^{) sn( Xl ) ^ sn{ Xl ) ^ 2sn(x 1 )J V ' 

We will consider m — 1 > c 3 > 0. For these values of m, the potential in H m is confining, i.e. it diverges 
to infinity in the boundaries. Energy levels of H m are doubly degenerate including the non-zero ground 
state energy £ m = (m + — c i- Doublet of ground states is formed by and where explicitly 



ip + = y / dn(xi)sn(xi)2 +m 



\ — 

1 — cn(xi) \ 2 
1 + cn(xi) 



( 



1 / 1— cn(xi) 
n — C 3 + m \ t, 7 f 

2 ° \l+cn(xi) 



1 / 1— cnfaii) 

V -\/-2- c 3- m {T+^T) 



(56) 



The mass M(xi) = dn (xi, 1) = sech xi coincides with the mass term used in [42] . 
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The vector complies with the required boundary conditions. The action of j m on (|56p preserves regularity 
of wave function provided that m is integer (semi-integer) and C3 is semi- integer (integer) . The supersym- 
metry ()42|) of H m is spontaneously broken; the supercharges (|53p and ia 3 Q^ do not annihilate ground 
states, but interchange them mutually. 

The superpartner Hamiltonians H m and H m+ i are spectrally almost identical up to the energy level 
£ m which is missing in the spectrum of H m+ \. Hence, the extended operator % m in (|44p has four- fold 
degeneracy of the energy levels up to the ground state which is doubly degenerate. The action of the 
operators Q$ keeps invariant the subspaces with fixed m. The operators jffl annihilate the ground 
states of 7i m , see in Fig. 2 for illustration. 



in+4 - • o ■ □ 




m 771 + 1 



Figure 2: (Color on-line) Left: the spectrum of T-L m . The operator Q m (solid line) keeps the doublets 
of H m (dots) or H m+ i (squares) invariant. J7m transforms eigenstates of H m into eigenvectors of H m+ \ 
and vice versa. Right: probability density of the ground state (blue dashed line) and the first excited 
states (red dashed line) of 1+ ^ 3 H m . Solid black line corresponds to M(x\) = dn(xi) 2 , dashed black line 
corresponds to the logarithm of the potential. We fixed m = 2, C3 = \ and k = 0.9. 



5 Discussion 
5.1 Dirac oscillator 

The quantum settings presented in section [3] were based on the potential algebra so(3). Let us make few 
comments on the quantum systems whose integrals of motion form either the oscillator algebra or so(2, 1). 
When ci = 0, one of the solutions of ([5]) corresponds to g\ = 1, gi = and 33 = C3X1. Substituting these 
coefficients into ([T]), we get the Hamiltonian 

h D = ia\d Xl + <7 2 h C302Z1, x\ > 0. (57) 

which coincides with the radial part of the well known Dirac oscillator Hamiltonian introduced by Moshin- 
ski and Szczepaniak [45 1, It has integrals of motion in the following form, 

J± ^^_(^ ±C3I1± -). (58) 

The eigenfunctions of ()57p are required to be vanishing at x\ = 0. The system has bound states as long 
as C3 / 0. They can be generated from the lowest weight vectors 




i;-=e mX2 e^x^ m [-iJc 3 (-^ + m),V2c 3 x 1 } , (59) 



13 



and comply with the required boundary condition provided that C3 < and m is a semi-integer lower or 
equal to |. 

In general, the solution of the stationary equation with (|57|) are not determined uniquely by the 
requirement of square integrability. The singularity of the potential at x\ = causes that the square 
integr ability is not sufficient to determine the wave functions. Our requirement of regularity of the wave 
functions at the origin removes the ambiguity by fixing one of the self-adjoint extensions of hp- In [49J, 
different way of dealing with singular potentials was discussed. The potential of Calogero Hamiltonian 
was regularized by the complex shift of the coordinate. It violated its hermicity, however, the complex 
shift of coordinate allowed for additional set of solutions and spectrum remained real (see also [50] for 
a recent reference where the technique was used). This approach might be applicable in the context of 
Dirac Hamiltonian (|57[) as well as in other quantum systems described by ([1]) with singular coefficients 
51, g 2 and g 3 . 

5.2 so(2, 1) and BTZ black hole -like system 

Finally, let us discuss the case where the quantum system possesses integrals of motion that form so(2, 1). 
We shall consider x 2 as the time-like coordinate. In order to do so, we perform the substitution x 2 — > —it. 
It makes the replacement dx\ — > —dt 2 in the line element ds 2 = g lLV dx il dx v and, hence, the sign of the 
lower component of g^ v gets changed accordingly. This substitution transforms Dirac Hamiltonian (|70p 
into (fTTj) (the flat space metric is now rj^ u = diag(l,— 1)) while keeping the condition (fT3|) unchanged. 
As an example, let us consider the space with the metric tensor in the following form 

9ilv = diag (N- 2 , -N 2 ) , N = N( Xl ), (60) 

where we leave the function iV unspecified at the moment. Let us comment that metric tensors of similar 
type appear in description of black holes [51J. The Dirac Hamiltonian (|71|) associated with (|60|) reads 

fiD = iNa%d Xl + + icri^- + c 3 a 2 N' . (61) 

It has integrals of motion (|12p , 

J± = fe* (Nd Xl =f N'd t T I T ||) ■ (62) 

They, together with J3 = dt, form the so(2, 1) algebra provided that the components of the metric tensor 
([60]) solve (fl~3|) . We find very interesting that the solution can be found where the function N acquires 
the following form 

iV 2 = _l +x 2_ ( 63 ) 

The metric (j60|) then coincides with the one of BTZ black hole [52] as long as the black hole has unit 
mass and zero angular momentum. The particle lives on the straight line with fixed angular coordinate 
which terminates in the singularity. This coincidence suggests that appropriate (higher-dimensional) 
modification of our current algebraic approach could be particularly helpful in the analysis of Dirac 
fermions in the gravitational background of BTZ black hole. 

There appeared remarkable results recently where the geometrical properties of black-holes were en- 
coded into the spatial part of the metric tensor. This observation enabled to consider black hole simula- 
tions in the strained graphene [53] , |54| . Analysis of these systems in the context of the potential algebras 
could be very interesting. 

5.3 Outlook 

In the current work, we considered potential algebras of the (2 + 1)- dimensional Dirac operator. The 
explicit form the operator was found that possessed Lie algebra so(3), so(2, 1) or oscillator algebra of 
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integrals of motion. We showed that the algebraic framework can be very useful in the analysis of various 
quantum systems. 

The structure of the potential algebra, its rank in particular, was fixed by the ansatz and the 
explicit form of the Dirac Hamiltonian (pQ). It would be possible to extend the current approach to 
higher-dimensional settings by considering potential algebras of higher rank. Such algebraic structures 
could be used effectively in the analysis of quantum systems that live in spaces with non-trivial geometry. 
It might be applied in construction of the new shape-invariant two-dimensional Dirac Hamiltonians. As 
it was suggested above, it could be also useful in the analysis of Dirac particle in the gravitational 
background of BTZ black hole. 

In the considered framework, the operators J± were required to commute with the Hamiltonian, i.e. 
they reflected degeneracy of its spectrum. The concept of dynamical algebras (or spectrum generating 
algebras) represents a natural extension of the current approach. Relaxing the condition [J±, hp] = 0, we 
could consider a broader class of quantum systems. A possible hint on how to modify the commutator is 
provided by the shape-invariance in the nonrelativistic quantum mechanics, where -£f~ +1 g m — QmH^ = 
r(m)q m . Having in mind the reduction ([32]) and (f33|) . the commutator is suggested to be [ho,J±\ ~ J±. 
The ansatz for the operators J± should be extended appropriately to the matrix operator with generally 
non-vanishing antidiagonal components. 

It is worth mentioning in this context that a different implementation of dynamical symmetries to 
nonrelativistic systems was introduced in [55J. Dynamical algebra of Poschl- Teller system was treated 
as a potential algebra of a suitably modified Hamiltonian. Implementation of this approach to Dirac 
operators could be fruitful as well. 



Appendix A 

We shall show that for the representations of so(3) generated by (HJ) from the lowest weight vectors, the 
coupling constant C3 has to be integer (half-integer) as long as m is half-integer (integer). Only in this 
case, the representations can be finite- dimensional. 

To show this, let us suppose that <?2 > and fix the function g\ in the following way 



9i\Jg\ - 1 
91 = — — • 

Then the lowest and highest weight vectors can be computed explicitly as 



(-l + gl)± \^g^w^J v V 9 V 



where the actual value of the constants and (3^ is not important at the moment. Next we define a 
matrix function G, 

92 \ g 2 2 w 1/z ) 

We use this matrix to transform the ladder operator J + into the more convenient form, 

1 n J3 



G- L J+G = ie^g 2 ^-l + gl[^d Xl +^) (64) 

\9 2 92. 



The function ij) can be rewritten correspondingly as 



m-l+2c 3 „„2c 3 -l 
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Let us compute how does J + act on the generic wave function £ = e mX2 G \^ 2W J. Using ([MP , we can 
calculate immediately 



J+e inx 2 Q f 92 r ' wS =ie i(n+l) G f -(n + r-s)w s 1 g r 2 1 + g r 2 +1 w s (n + r) ^ 



(66) 



This formula helps to understand qualitatively the structure of the (J+) fe £ (with k being a positive integer); 
we can find coefficients at the highest and the lowest power of g 2 - The coefficient of the term g 2 ~ k w s ~ k 
is (-l) fe Ui=o( n + l-l + r-s) while the coefficient of the term g r 2 +k w s is n?=ci( n + r + 2 0- 

Now, let us identify the nonzero element of £ with the upper component of (|65p. i.e. we fix n = m, 
r = m — 1 + 2c3 and s = 2c^ — 1. The requirement that is annihilated by a specific power of J + implies 
that the above derived coefficients have to vanish for a specific value of k. This happens provided that 

2m - 1 £ Z < 0, -(m + c 3 ) + ^ G Z > 0. 

The first relation restricts m to be integer or half-integer while the second relation tells that C3 has to be 
half-integer or integer, respectively. 



Appendix B 

Let us review briefly how the Dirac Hamiltonian for mass-less particle in the curved space can be con- 
structed. We will consider two distinct scenarios that are distinguished by the sign e in the metric tensor 

_ __ ( 5ii (a?i) 
9 ^~\ eg 22 (x l ) 

We suppose that gu{x\) > and g 22 (xi) > 0. For e = — 1, we deal with (1 + 1) dimensional space-time, 
considering the coordinate x 2 to be time-like. When e = 1, both coordinates are space- like and the metric 
describes curved euclidean space. In both cases, Dirac equation can be written in the following form [23], 

*/• (0„ + t/> = A^, = 1,2, (67) 

where A corresponds to the mass of the particle for e = — 1 whereas it represents eigenvalues of Dirac 
Hamiltonian for e = 1. In the latter case, (|67p can be identified with the stationary equation. is spin 
connection and 7^ are generalized gamma matrices which satisfy {7^,7^} = 2g fiU . These quantities can 
be computed with the use of tetrad formalism. The tetrads are related to the metric tensor in the 
following manner, 

e x ll e x l/ r lab = 9iMUi 

where the quantities e" and e% relate locally flat coordinates y a with the curvilinear coordinates at 
a given point X, e Xfi = (d Xll y a ) (X) and = {d ya x^) (X). There holds e^e^ = 5 b a . Flat-space metric 
is identified as 77^ = diag(l,e). The matrices 7^ are defined as 

7 ^ = e£ 7 a where { 7 a , 7 6 } = 2r] ab . (68) 

The spin connection is defined as in terms of tetrads and affine connection is T^ u as follows 

1 

t 
where 

< = e a Xx g Xr (d^e b XT - T^) , r*, = \g x ° (d v9afl + d^ au - d a9flv ) . (69) 
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In our present case, the tetrads can be fixed as 



/9n \ x M _ / Tfrr 

^ I ' " I 



922 



Here, the upper index denotes row and the lower index the column of the matrix. Non-zero elements of 
the affine connection are = T\ 2 = s-^, rf 2 = = and 



, ,12 _ , ,21 _ 922 
U> 2 — —U! 2 



2 a/511522 

The matrices 7 a are defined as 71 = o"i, 72 = i~2~cJ2 and satisfy {7a, 7&} = 2r]ab- The generalized gamma 
matrices are then 

7 s1 = \/5lT 1 cri, 7 X2 = £^\/522~W 

Finally, the spin connection is 

* (^) 922 

Q Xl = 0, U X2 = -o- 3 . 

Dirac Hamiltonian for e = 1 acquires the form 



&r = ^1 -7=^ + 7-^7= )+^^=, (70) 

/5IT 45r 2 2A/5TT/ a/522 



while it has the following explicit form for e = — 1 



/4T 1 = "1 ( ~^d Xl + 7^^) - a 2 ^. (71) 

4022^511/ V^22 

We transform the operator /if^ 1 to the form which is manifestly hermitian with respect to the standard 
scalar product in the following manner 

h D = (detg^H D (detg^)- l 4 = ^d xi + -^d X2 - ^± (72) 

y/9u V92^ 4^f 
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